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ABSTRACT 

FUor outbursts in young stellar objects (YSOs) are the most dramatic events 
among episodic accretion phenomena. The origin of these bursts is not clear: 
disk instabilities and/or disk perturbations by an external body being the most 
viable hypotheses. Here, we report our VLT/SINFONI high angular resolu¬ 
tion AO-assisted observations of 2MASS J06593158-0405277, which is under¬ 
going a recently discovered FUor outburst. Our observations reveal the pres¬ 
ence of an extended disc-like structure around the FUor, a very low-mass com¬ 
panion (2MASS J06593158-0405277B) at ~100 au in projection, and, possibly, 
a third closer companion at ~llau. These sources appear to be young, dis¬ 
playing accretion signatures. Assuming the components are physically linked, 
2MASS J06593158-0405277 would then be one of the very few triple systems ob¬ 
served in FUors. 

Subject headings: binaries: close — circumstellar matter — stars: formation — stars: 
pre-main sequence — stars: individual (2MASS J06593158-0405277) — techniques: 
high angular resolution 
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1. INTRODUCTION 

FU Orionis-like outbursts (hereafter FUors) are the most dramatic events among 
episodic accretion phenomena in YSOs, displaying a mass accretion rate (M acc ) increase 
by several orders of magnitude (up to M acc ~10” 4 M Q yr _1 ; see e.g. Hartmann & Kenyon 
1996; Audard et ah 2014). As a result, the YSO brightness is boosted by 5 magnitudes 
or more in the optical (e.g., Herbig 1977). These bursts have a dramatic impact on the 
YSO circumstellar environment (Baraffe et al. 2012; Vorobyov et al. 2013; Audard et ah 
2014). The rise times of these outbursts are usually short (from several months to a few 
years), whereas the decay timescales range from years to several decades. Spectra of FUor 
progenitors have only been obtained in two cases, both indicating that they are low-mass 
YSOs (e.g., Herbig 1977; Miller et ah 2011). 

The origin of bursts in FUors is not well understood, and disk instabilities, disk 
fragmentation, and/or disk perturbation by an external body are the most viable 
hypotheses (see review by Audard et ah 2014, and references therein). Several models 
with thermal or gravitational plus magneto-rotational instabilities are able to mimic the 
main characteristics of the outbursts (see, e.g., Clarke et ah 2005; Zhu et ah 2010a, b). 
Notably, the recent version of the disk gravitational instability and fragmentation model by 
Vorobyov & Basu (2015) reproduces the physical properties including the episodic nature 
of the bursts. Aside from internal instabilities, external triggers, namely close encounters in 
binary or multiple systems, can produce FUor bursts (Bonnell & Bastien 1992; Reipurth & 
Aspin 2004; Reipurth et ah 2014). 

FUors are rare transient phenomena, and they have usually been detected and studied 
in detail several months or even years after the beginning of the outburst. 2MASS 
J06593158-0405277 is a YSO undergoing an FU Ori-type outburst (Maehara et ah 2014; 
Hillenbrand 2014). The outburst was detected on the 23rd of November 2014 (Maehara 
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et al. 2014), and previous observations indicate that the object began its slow rise at the 
end of 2013 and reached its maximum at the end of 2014 (Hackstein et al. 2014; Kospal 
et al. 2015). From archival photometric data, Kospal et al. (2015) identify the progenitor 
as a low-mass T Tauri star (0.75 M 0 ) with a circumstcllar disk of 0.01-0.06 M e , and an age 
of ~ 6xl0 5 yr. 


2. OBSERVATIONS AND DATA REDUCTION 

2MASS J06593158-0405277 was observed with ESO/VLT UT4 on 2015 January 25 
with the near-infrared integral field spectrograph SINFONI (Eisenhauer et al. 2003) in the 
J (1.11.4 /mi) and K (1.95-2.45/mi) bands with a spectral resolution of R ~ 2000 and 
4000, respectively. The observations were performed with the highest spatial sampling (12.5 
milliarcseconds - mas - pixel scale, and field of view - FoV - of O'/S x O'/S) and in AO-assisted 
mode (using the target as a natural guide star) with a mean optical seeing of ~ O'/S, which 
provided a full-width at half maximum (measured on the point spread function - PSF) 
of ~ 50 mas and ~ 70 mas in the J and K band, respectively. The detector integration 
time (DIT), number of sub-exposures per frame (NDITS) and exposures were 30 s, 4, and 
5 (J band), and 15s, 5, and 7 (K band), for a total integration time on-source of 600s 
and 525 s in the J and K band, respectively. An additional photometric standard star 
(HD 42655) was observed for flux calibration of the data cubes and removal of telluric 
features from the spectra. Data were reduced using the SINFONI data reduction pipeline in 
GASGANO (Modigliani et al. 2007), that is, dark and bad pixel masks, flat-held corrections, 
optical distortion correction, and wavelength calibration with arc lamps. 

Two spectra were extracted from each data cube using a 5x5 pixel area centered at 
the position of the peak flux of the FUor and its companion (B, see Sect. 3). Due to the 
FUor brightness and its spatially extended disk-like emission (see Sect. 3), a background 
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spectrum was extracted close to the companion position and subtracted from its spectrum 
to remove the contamination from the FUor source. 

Additionally, continuum images (at 1.25 and 2.2 pm) and continuum-subtracted Pa/3 
and Bry images were created following the same procedures as in Beck et al. (2008) and 
Garcia Lopez et al. (2013). 

As no PSF standard was observed, to remove the FUor’s PSF from the continuum and 
evaluate the asymmetry of the bright central core (see Sect. 3), we rotate the continuum 
images by 180°, match the PSF and subtract the rotated images from the original ones in 
both J and K bands. 


3. RESULTS 

As illustrated in Fig. 1, SINFONI data cubes in both bands show three main features 
around 2MASS J06593158-0405277: i) a faint companion (hereafter 2MASS J06593158- 
0405277B or B ) located at 227 mas (or 102 au in projection at a distance of 450 pc 1 , 
see Kospal et al. 2015), southeast from the central source (see left panels of Fig. 1); ii) 
an extended and elliptical disk-like structure around the FUor, (size ~0'/2 or 90 au at a 
distance of 450pc), possibly scattered light from the disk or circumstellar nebulosity (e.g. 
Goodrich 1987, see also left panels of Fig. 1); in) a bright asymmetric central core of the 
FUor source, which has a ‘bump’ towards the southwest, clearly detected after subtracting 
the FUor’s rotated images (see central panels of Fig. 1). Moreover, both Pa/3 and Bry 
lines in emission are detected at the position of 2MASS J06593158-0405277B and near the 
continuum emission of the ‘bump’, as shown in the right panels of Fig. 1. 

1 The distance to the object is highly uncertain as well as its association with the molecular 
cloud L 1650, d~2.3Kpc (Kim et al. 2004) 
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Fig. 1.— Left: SINFONI J- (upper panel) and A"-band (lower panel) images of 
2MASS J06593158-0405277 showing the disk-like feature, the “bump” position, and the com¬ 
panion (B). Center: Difference J- (upper panel) and A"-band (lower panel) images showing 
the ‘bump’. The PSF was removed by subtracting the same image rotated by 180°. Right: 
Pa/3 (upper panel) and Bry (lower panel) continuum-subtracted gray-scale images. The same 
continuum contours are shown as in the left panels. 
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3.1. The Companion: 2MASS J06593158-0405277B 

On the basis of Gaussian fitting of both J and A"-band data cubes, we infer the 
following coordinates for B\ R.A.(J2000)=06:59:31.599, Dec.(J2000)=-04:05:28.01, that 
is, its separation from the FUor is Aa=0'.'15, Ah=-0 / .'17, with a P.A.=131?4. Our PSF 
fitting allows us to derive J and A'-band magnitudes of 14.8±0.4mag and 12.9±0.2mag, 
respectively. The photometric uncertainties mostly arise from the bright background due 
to the scattered light from the FUor. In the upper panel of Figure 2, we show the flux 
calibrated spectrum in the J and K bands. The spectrum was obtained after subtracting 
the FUor scattered light in the background (see Sect. 2). The spectrum shows several 
broad molecular absorption bands from H 2 0 (1.1, 1.4, 1.9, and 2.4/im), CO (bandheads 
v=2-0, 3-1, 4-2, 5-3 from 2.29pm in the K band), and, possibly, FeH in the J band, 
as well as atomic absorption lines (Ca I and Na I) in the K band. These absorption 
features have a photospheric origin, and are distinctive of late M-type stars (see e.g., 
Cushing et al. 2005; Rayner et al. 2009; Scholz et al. 2009). We use near-IR spectra of very 
low-mass and brown-dwarfs YSOs from Scholz et al. (2012) to constrain the spectral type 
of B. The inferred spectral type (SpT) is between M5 and M9, that is SpT=M7±2, with 
T e ~2900±300K (Luhman et al. 2003). To get a rough estimate of the stellar parameters, we 
assume that B is coeval to the FUor, thus with an age between 5 x10 s and 10 6 yr (Kospal 
et al. 2015). We adopt the evolutionary tracks from Baraffe et al. (2015), inferring 
M* = 0.1 ±o;o 8 M q and A* = l.Oi^gR©. 

We use the J — K color excess of 2MASS J06593158-0405277B and an SpT=M7 (with 
(J — AT) 0 =1; see e.g., Scholz et al. 2012) to estimate both upper and lower values of the 
visual extinction (A v ) towards the source, adopting the extinction law from Rieke & 
Lebofsky (1985). Indeed, by assuming that the observed E( J — AT) is only due to the visual 
extinction, we get an upper limit of Ay=5 mag from Ay = [( J — K ) — ( J — A") 0 ]/0.17. This 


is only an upper limit, because this source is young and thus part of its IR excess originates 
from veiling (r A = F ex ^/F*\). The IR excess is thus: 


E(J - K ) = (J - K) - (J - K )o + 2.5 log ( (1) 

\l + r K J 

We cannot derive both rj and tk values. However, by assuming rj=0 and fitting 
an SpT=M7 template, we then get r^~2.6, and we can infer a lower limit for the visual 
extinction of Hy=2.4mag. 

The spectrum of 2MASS J06593158-0405277B (see upper panel of Fig. 2) also displays 
H I lines in emission, namely the Pa/3 line at 1.282/rm (F = (2.4±0.1)xl0 _15 ergs _1 cnr 2 ) 
and the Bry line at 2.166/im (A = (6±2)xlO” 15 ergs -1 cnr 2 ). These lines are associated 
with ongoing accretion (see, e.g., Muzerolle et al. 1998; Natta et al. 2004; Caratti o Garatti 
et al. 2012), and can be employed to derive an estimate of the mass accretion rate (see 
e.g., Natta et al. 2004, 2006). With the empirical relations by Muzerolle et al. (1998), this 
translates into an average accretion luminosity of L acc =0.005-0.01 L Q . The mass accretion 
rate can then be expressed as M acc = 1.25 x L acc U t /GM* (see e.g., Caratti o Garatti et al. 

2012) . Inferred values are between 4xl0 -9 and 8xlO~ 10 M 0 yr _1 , consistent with accretion 
rates found in young brown dwarfs (Natta et al. 2004; Mohanty et al. 2005; Joergens et al. 

2013) . 


3.2. The Disk, the Bump, and the FUor 

J and K band continuum images in Fig. 1 (left panels) display an elongated disk-like 
structure extending up to ~200 mas (or ~90 au) from the central source. Elliptical isophotes 
were fitted to the disk-like feature using the iraf STSDAS task isophote, starting from radii 
larger than the central PSF (i.e. ~50mas and ~70mas in the J- and A'-band, respectively). 



F a (erg cm 
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Fig. 2.— SINFONI J- and A"-band flux calibrated spectra of 2MASS J06593158-0405277B 
(red solid line, upper panels) and of the FUor (black solid line, lower panels). The most 
prominent features are also labeled. 
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The isophotal fit in the J-band provides us with a semi-major axis P.A. of 52°±3°, and a 
ratio of the minor to major axes of 0.88±0.05. Assuming circular symmetry, the system axis 
inclination with respect to the line-of-sight is i=28°±6°. Results of the K band analysis are 
less accurate but consistent with the previous ones, yielding P.A.=48°±10° and i=22°±12°. 
Our analysis thus indicates that this disk-like structure is almost face-on. 

Another notable circumstellar feature around 2MASS J06593158-0405277 is the ‘bump’ 
detected in continuum and line emission images (see Fig. 1, central and right panels). This 
feature is not spatially resolved in the J and K band images, but it is clearly detected 
after subtracting the rotated images (Fig. 1, central panels) as an arc-shaped elongated 
feature, peaking at P.A~233° and ~225° in the J and K bands, respectively, and with 
mj~12.9mag and m^~10.6 mag. In the continuum-subtracted Pa f3 and Bry images 
(Fig. 1, right panels), the ‘bump’ emission appears as a point-like source, well fitted with a 
Gaussian profile. Its coordinates are R.A.(J2000)=06:59:31.588, Dec.(J2000)=-04:05:27.86, 
i.e. it is at Aa=-15mas, Ah=-20mas from the FUor, that is Ar=25mas (or 11.25au), 
with a P.A.=233?1. A rough estimate of the Pa (3 and Bry line fluxes gives ~7xlCR 15 
and 2xlCR 15 ergs" 1 cm" 2 , respectively. The detection in both continuum-subtracted line 
images, the fact that these features are fitted with Gaussian profiles, and the presence of 
continuum emission in the PSF-subtracted J and K band images strongly suggest that the 
‘bump’ is an additional closer young companion, which might have triggered the burst (see 
e.g., Reipurth & Aspin 2004; Reipurth et al. 2014). Moreover, in both J and K bands, the 
‘bump’ appears brighter than 2MASS J06593158-0405277B, hinting that, if it is a YSO, it 
might be more massive than the latter. However, as the ‘bump’ is not spatially resolved in 
the continuum, other scenarios cannot be completely ruled out for the continuum emission, 
such as an uneven disk illumination (see, e.g., Dullemond & Monnier 2010). 

Finally, the lower panels of Fig. 2 show the FUor flux calibrated spectrum 
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(mj=8.63±0.05mag and mif=7.45±0.04mag), which only exhibits molecular and atomic 
features in absorption, as usually observed in FUor spectra, namely H 2 0, CO, H I (Pa (3 and 
Bry), Na I, Ca II, K I, Mg I, and Fe I. From these features we infer a spectral type between 
K9 and Ml, typical of FUor NIR spectra (see, e.g., Audard et al. 2014). Notably, this 
spectrum closely resembles that of HBC 722, a well studied FUor (see, Miller et al. 2011). 


4. DISCUSSION 

The detection of binaries or higher order systems in FUors is potentially of great 
importance, as interaction with a companion could be a triggering mechanism, although, in 
principle, there might be more than one. For example, a couple of FUors (V1057 Cyg and 
V1515Cyg) are probably single systems (see Audard et al. 2014, and references therein), 
thus, in this case, internal instabilities seem to be the most straightforward process to 
consider. On the other hand, a large number of FUors and FUor candidates are binaries, 
and Bonnell & Bastion (1992) first proposed that FUor events may be triggered by the 
passage of a close companion. Later on, Reipurth & Aspin (2004) suggested that these 
objects might be newborn binaries, which arise from the breaking of unstable triple 
systems (see also Reipurth et al. 2014). As a consequence of the dynamical interaction, one 
object is ejected, typically the lowest mass component, whereas the remaining two objects 
spiral in toward each other perturbing their disks. The bursts would be then triggered by 
the close passage (~10 au) of the nearest companion. 

This scheme is a plausible explanation for what we are observing in 2MASS J06593158- 
0405277. The central object is a FUor, whose burst might have been triggered by a close 
companion (the ‘bump’) located at ~llau, and there is a more distant (~100au) and 
less massive object, which could be the third (ejected?) companion of a non-hierarchical 
system (e.g., Reipurth et al. 2014). Notably, the present configuration is precisely what 
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would be expected from the dynamical evolution of a newborn triple system, where the least 
massive of three bodies is kept away from the two other bodies except for brief periastron 
passages, which eventually lead to the system dissolution (see, e.g., Reipurth & Clarke 2001; 
Reipurth et al. 2010). At an age of ~1 Myr the observed system is statistically unlikely to 
be stable, and it will eventually eject the low-mass component after a periastron passage. 
This could be the first time that such scenario has been observed. Indeed companions in 
other FUors have been detected (e.g., FU Ori, Z CMa, L1551IRS5, RNO 1B/C, AR 6A/B; 
see e.g., Reipurth & Aspin 2004; Audard et al. 2014), but their relatively wide separation 
(from several tens up to hundreds of au) suggests that these sources might be outlying 
companions. 

In the case of 2MASS J06593158-0405277, the fact that both the ‘bump’ and source 
‘B’ are young objects excludes the possibility that they might be background or foreground 
stars. Instead it is much more likely that they belong to the same star forming region, 
which is quite small and isolated, being composed of a few TTauri stars (Kospal et al. 
2011). Their proximity to the FUor and the fact that the observed radial velocities of the HI 
lines are very similar in the three objects (average u r .~7±7kms~ 1 with respect to the local 
standard of rest -LSR- and not corrected for the molecular cloud velocity) suggests that 
these sources might belong to the same system. Provided that the objects are physically 
linked, 2MASS J06593158-0405277 would then be one of the very few triple systems ever 
observed in FUors. That said, further high-angular resolution observations are required to 
confirm the nature of the ‘bump’ and the multiplicity of the 2MASS J06593158-0405277 
system. 
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